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ABSTRACT 



The propagation of sound in a laboratory model 
consisting of a surface layer of water overlying a thick 
bottom layer of sand was experimentally investigated. The 
variations of pressure amplitude as a function of receiver 
depth, range, and transverse location were measured for the 
source at constant depth in the surface layer of water. 
Three frequencies for which only the lowest mode propagates 
were used in this investigation. The measurements of pres- 
sure amplitude as a function of depth and range were found 
in good agreement with normal mode theory, when absorption 
in the bottom was taken into consideration. 
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I. PROPAGATION OF SOUND IN THE HOMOGENOUS LAYERED MODEL 



The transmission of sound in the ocean is dependent on 
several environmental factors. The most important of these 
include the depth, the nature of the bottom, sound velocity 
structure within the ocean, and the shape of ocean bottom 
and surface. The purpose of this experiment was to investi- 
gate the propagation of sound in a two layered laboratory 
model, composed of sand and water. In this study the bottom 
and surface of the model were made flat and experimental 
results were compared with a simple theoretical model. 
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II. PROPAGATION OF SOUND IN SHALLOW WATER 



A. INTRODUCTION 

The propagation of sound in shallow water is a compli- 
cated problem, involving multiple reflections of the sound 
between the surface and the bottom, so that there is consid- 
erable phase interference among the various multipaths. At 
relatively long ranges, the sound field is mainly due to the 
multiple reflections. The transmission of sound in the 
shallow water also depends on several environmental factors: 
nature and slope of the bottom, depth of water relative to 
wave length, roughness of surface and bottom and temperature 
gradient. As a consequence of the multiplicity of reflec- 
tions from the surface and the bottom and the environmental 
factors, the predictions of sound transmission in shallow 
water can be much more complicated than for deep water. Ray 
theory can be used for the sound field at short range when 
the depth of the water is relatively large. But, ray theory 
has limited applications. It is inapplicable in a shadow 
zone. Ray theory also cannot be used for low frequencies 
when the wavelength of sound becomes comparable with the 
vertical scale of sound velocity variation. Therefore, in 
practice one ultimately has to resort to a wave solution of 
the problem. 

B. NORMAL MODES IN TWO HOMOGENOUS LAYERS 

Normal mode theory provides an exact solution to the 
wave equation and is also useful for dealing with depth 
dependent sound profile particularly those leading to 
the trapping of sound. There are various approaches to 
the normal mode theory but we will use the 
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eigenfunction-expansion technique, which is relatively 
simple. In this presentation we will follow mainly [Ref. 1] 
and [ Ref . 2 ] . 

Consider the case of two homogeneous layers with a 
pressure release surface above as shown in Figure 2.1 . 



Pressure release 




Figure 2.1 Two homogeneous layers with pressure release 
surface . 

Let and f , C 2 be the constant densities and sound 
speeds in layer 1 and 2 respectively. Assume there is a 
point source at z = d in layer 1. The wave equation for a 
homogeneous medium can be written as • 
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V p - t a 2 p _ s 

~ Par? ( 2 -d 

where (Tis a source function and is defined by 

cr = 6Q 

a t 



Where Q is source strength 

In our solution, we assume that a point source radiating 
a single frequency can be described mathematically by 



cr 



4TTA S( r - 




t 



(2.2) 



Where S ( r - r. ) is the Dirac delta function, r is the 
position vector. A is the amplitude factor which is 
assumed to be unity, and co is the angular frequency of 
the source. Solution to the wave equation is of the form 



PC r t ) = Pd-je'"' 

(2.3) 

Substituting equation (2.2) and equation (2.3) into 
equation (2.1) we get the helmholtz equation 



Vp 



+ 



2 

k P 



- 41T Sir - r ) 



(2.4) 



13 



where k is wave number defined by 



k = 



Go 

C 



(2.5) 



The delta function for the point source at r = d is 



S( r-r ) = _J_ 5(r)6(z -d) 

2TT r 



( 2 . 6 ) 



In cylindrical co-ordinates equation (2.4) can be written as 




2_<5(r)<$(z-d)-e yw/ 

r (2.7) 



This is an inhomogeneous equation. First, let us 

consider the homogeneous equation. 




A-(rd_) 
3r ar / 



ar 



+ 




k* ]f = 0 



( 2 . 8 ) 



and use the separation of variables technique. 

We assume that the sound energy is trapped in the upper 
layer and solution to equation (2.8) can be represented by 
a set eigenfunctions. 



fov/t) 



e' V T RO-)Z(z) 

n n 



(2.9) 
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Substitute equation (2.9) into equation (2.8) to get 




or 




= 0 (2.10) 



f J_J_ d_/ 


'rdR,\l f 


1 r d r ' 


■ dT 2J + 1 



2 

dZ „ 



Z„d 



+ 



] = 



0 



( 2 . 11 ) 



This is only possible if each part in the square 
brackets in equation (2.11) is equal to a constant hence it 
can be written as 



JLii-rrd&A = 
FL r dr^ dr/ 




( 2 . 12 ) 



-L + 



= k. 



(2.13) 



where k n is separation constant. 

The boundary conditions on the depth functions are 



Z,„(0) = 0 



(2.14) 
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(2.15) 



Z (H) = 


Z (H) 


in 




if dZ,„\ 


1 f dZi/i\ 


P K dz > H 


p k d Z 'H 



lirn Ijz) = 0 

Z-+OO 



Rewrite equation (2.13) in the form 



or 



where 



2 2 
d Z n , f k 






+ 



( 



K^n 



= 0 



4. k Z — 0 

d 1? zn n 



k = (co f k 2 
zn {-cl) ~ 



In order to satisfy the boundary condition 
goes to zero when Z — ► oo , the solution Z tn 
exponentially with depth. The only way that Z, n 



(2.16) 



(2.17) 



(2.18) 



(2.19) 

the pressure 
must decay 
can satisfy 
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the boundary conditions at the surface at Z = 
be sinusoidal function of depth. Therefore 



His to 



U 

II 


A S i n 
n 


k z 

zn 


0 H 


(2.20) 


Z = 

zn 


.8 z 

r-s, “n 

B*e 

n 




Z< oO 


(2.21) 


2 


2 


2 






R, 


ii 


-(f) 




(2.22) 


Waveguide propagation 


requires 


that k n 


and must 


be real 



hence 



< k 



< 
n ~~ 



co 

C, 



(2.23) 



Application of boundary Conditions at Z = H results in 



A - sin k- h 

zn 




(2.24) 



A 



r\ 



1 k cos k H 




zn 




(2.25) 
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which gives 



B„ = A/ W sin k z h 



(2.26) 



The depth function, therefore, can now be written as 



Z Cz) = A sin k • z 

i n n zn 



0 < Z < fi (2.27) 



Zfz) = 

in 



A sin k He 

n zn 



-Pjz-H) 



2 < « (2.28) 



These depth functions have to be orthogonal. Since there 
is a discontinuity in the slope at the boundary Z = H, we 
need separate orthogonality constant for each layer. We 
assume that the constant is unity for the layer 1 and 7S for 
the layer 2. 



Z„(z) = I 






7Z !n (z) 



0< Z< H 



H < Z < co 



(2.29) 



In order to evaluate "3 multiply equation (2.13) by Z 



2 2 



Z m dZ , (k-k)zz _ 

^ ^ n n m 



0 



(2.30) 
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interchange the subscript m and n in equation (2.30) and 
take the difference of the two equations to get 



2 .2 



i (Z dZ„ z dZ^ = fk-k m )ZZ n l2 31) 

Integrate both sides of equation (2.31) from 0 to °0 



f* d (Z dZn ZdZ m \dz = ( 
J dz'' m dz dz > 



2 2 , 
K - 



, oO 



(2.32) 



The orthogonality condition requires that, for m ^ n, 

the R.H.S of the equation (2.32) must be zero. Therefore 



f d-(z im ]!L«-Z w 4Lyz + 7 fuz, dZ,.z 2>) ^)dz = o 

J d z ^ dz dz' ' dz v dz dz ; 



(2.33) 



H 2 



[Z,_dZ,„ 


Z,„dZ,. 


1 + h ? Z„ d Z 


d Zj*l 
ZY> 


L d z 


d z 


J o 1 d Z 


d z J 



(2.34) 



The application of limits as well as boundary conditions 
gives 





0 (2.35) 
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This is only possible if 



i_7A = 0 
A 



or 



7 




A 

A 



(2.36) 



Therefore, equation (2.29) becomes 




0 < Z < H 



H Z < oo 



(2.37) 



We will normalize the depth eigenfunctions in order to 
evaluate A n , consider the m = n case 



Z 2 ( Z )d Z 
n 






(2.38) 



.H 2 

A Z Z dz 

/) i/j m + 



A 4 

Ja 



z Z dz 

2/7 2/7 



= 1 (2.39) 



V) - 



H Z 2 dz + A r°°z 2 dz 

in + / 2/7 



A 






(2.40) 
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Now 



r H 2 , H 

/ 2 - di - i 

0 v <L 



i_( 1 - cos 2k - z )dz 

2 zn 



(2.41) 



Integrating and substituting the limit we get 



r H 2 

/ Z dz - ir H L_ sin k • H -cos k • H 1 

1 ,n 2 L ~ 2n zn J ( 2 - 

0 /n 



42) 



and 



/ 



cr> 2 
2 » dz 



= rW^H^J-e^dz 



(2.43) 



or 



I 



oo 2 

Z. dz 



zn 



_l_Sin 2 k H 

« 



(2.44) 



Substituting equation (2.42) and equation (2.44) into 
equation (2.40) gives 



A. 



2 k 



zn 



kjd -coS(k^H)sin(k 7n H) -f/lfsinfkJH )tan(k„H ) 



\n J zn zn 

K (2.45) 
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now form a complete set of 



The depth eigenfunctions 
orthonormal functions. 

Mow we have to evaluate the Green' s function for the 
inhomogeneous wave equation 



Zz J_d_ 
" " r a r 




+ 




d 2 z 



0z 







2.5(r)5Cz-d) 

(2.46) 





2 

r 



S(r) S(z-d) 

(2.47) 



From equation (2.13), the second term in the square 
bracket is zero; therefore, equation (2.47) can be rewritten 
as 




_ JL r ) 8(Z -d ) 

r (2.48) 



Multiply both sides of equation (2.48) by Z w and 
integrate from 0 to co gives 




0 



<S(z-d)Z„d z 

(2.49) 



For m = n 



j d / r d B/A , 

r dr^ dr ' + 



k^R 

/) ^ n 



_ 2_S(r)Z (d ) 

r /) 



(2.50) 
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Solution to this homogeneous equation is 



( 2 ) 



R(r) = -j7TZ „(d)H (k r) 

'■-n J n 0 n 



(2.51) 



where for k^r >> 1 



( 2 ) , 

H (k„r) 






'ITk r 



The Hankel function H 0 (k^- r) describes the solution in 
the r direction. 

Therefore, the solution for the sound trapped in layer 1 
(0 < Z < H) is 



FKz.t) = -) 



-il 



27T A*sin(k-d)sin( k z ) + 4 ; 

L r n zn 

k„ r (2.52) 



Equation (2.52) is invariant under the interchange 
of Z and d and is valid within the region 0 < Z < H . 
The eigenvalues k 2n can be obtained by eliminating B n from 
equation (2.25) 



tan k H 

zn 



A 

A 6 



* ^ zn 



(2.53) 



Let 



X. - 



k-H 

zn 



b _ 



A 

A 



(2.54) 
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tan X 

n 




(2.55) 



tan X n 

x„ 



- 1 

bH V Co 2 ^ XT (2.56) 

c, 2 0/ H 2 



When C z > C, which is the case under consideration, 
equation (2.56) possesses real roots for k„ and X n , the 
later varying in the range (n - 1/2) 'TT < X ^ < n iT' . On the 
other hand, when C 2 < C, , the roots k n are complex numbers 
with negative imaginary parts. In the later case, there- 
fore, the factor exp(-jk >; r) in equation (2.52) implies 
horizontal attenuation. In the case of a fast bottom, there 
exist solutions which suffer no horizontal damping. This 
condition stems, of course, from the fact that for angles of 
incidence (measured from the normal) greater than critical 
angle, no power is transmitted into the bottom. Figures 
2.2, 2.3, and 2.4 shows the z dependence of the first three 
normal modes for 50 kHz, 120 kHz, 160 kHz respectively. 
Figure 2.5 shows the r dependence of the sound pressure. 
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Figure 2.2 Pressure amplitude vs depth for first mode. 
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Figure 2.3 Pressure amplitude vs depth for second mode. 
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Figure 2.4 Pressure amplitude vs depth for third mode. 
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Figure 2.5 Pressure amplitude as a function of range. 



C. CUT - OFF FREQUENCY 



Let us determine the cut-off frequencies of the normal 
modes. It follows from equation (2.56) that 



tan Xn 




-1 




bV "’ H UH)- x: 


(2.57) 



It is clear from the above equation that the nth mode 
cannot support propagation until 



2 

CO 




1 

c* 



?-*> 




(2.58) 



Therefore, the limiting frequency below which transmis- 
sion of the nth mode cannot take place is given by 



f = (2n - 1 )c, 

n ~ - v 

4H/ 1 - cf 
7 c? 



(2.59) 



where n = 1,2,3, 

This is an important relation between the layer thick- 
ness H and the cut-off frequency f n . It is shown in Figure 
2.6 for the first mode. 
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Figure 2.6 Cut-off frequency as function of layer thickness. 



D. PHASE SPEED 



The factor exp j (cot-k >; r+ tf/4 ) in equation (2.52) allows 
an obvious interpretation of k M , namely k n = c*yCpn where 
C pn is the phase speed of the nth mode. The phase speed is 
the speed that the wave front (surface of constant phase) 
propagates in r direction. The phase speed at the cut-off 
frequency has the value of the speed of sound in the 
bottom and decreases continuously with the increasing 
frequency towards the value of sound speed in the upper 
layer. From equation (2.19) we can write [Ref. 3] 



2 2 



k 



n 




(2.60) 



Multiply both sides by H gives 




(2.61) 



or 




or 




(2.62) 
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where 



oL 



n — 






Therefore equation (2.57) can be written as 



or 



ton X /? 

X„ 



-1/b 






tan Xn - 



- 8 



VcfV- x 



z 

n 



where 



6 = 

A c i 



From equation (2.65) 



'■Cif _ 1 



(t) 



(X 

n ~ 



6 



n 



(2.63) 



(2.64) 



(2.65) 



( 2 . 66 ) 



(2.67) 



( 2 . 68 ) 
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We know that C ^ n = k n H where k n = co/C pn . 
C pr1 is the phase speed of the nth mode. 
Therefore 



oo-H 



(2.69) 



From equation (2.62) we have 



co-H 



c V~* + X 

' n n 



(2.70) 



Hence 



C 

Pn 




(2.71) 



Substituting the value of from equation (2.68) we get 



C 

Po 




2 

6 



2 



1 + 6 2 cot 2 (X n ). 



(2.72) 



(n-l/2)1T <_ X ^ < n 1Y for the nth mode, and in this case 
k n and phase speed are real and propagation of the modes 
proceed without damping. X„ is a function of angular 
frequency and the frequency corresponding to any given 
value of X„ can be obtained. 
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y = cj-h 

27T C, 



H. 

X 



(2.73) 



= Cpn 

2 7T C , 

= X n Vi + g + fi-cot X, ’ 

27T6 

When X n is near its upper limit of nT(,then cot X — ► cO. 
Therefore is very large and from equation (2.72) 




Hence in the limit of very high frequency, phase speed 
approaches the sound speed . When X n is near its lower 
limits (n-l/2)TT then cot X — + 0 . 

2 i 

and C = cf 1+6 l 1 — C 

i [ J 2 

The phase speed, therefore, varies from the value C x at 
the critical frequency down to C, in the limit of very high 
frequency. Figure (2.7) shows the phase speed as a function 
of frequency for first three modes. 

E. GROUP SPEED 

In a medium in which the phase speed varies with 
frequency and if the pressure pulse from a point source is 
not a single frequency but covers a moderately broad 
spectrum, then the presence of the particular frequency 
component of the spectrum at a large range would occur at a 
time associated with the group speed (speed with which the 
energy associated with a particular frequency will travel). 
At frequencies co for which the phase is stationary, 
mutual interference will be negligible and these 
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be dominant 



frequencies 
of t and r 
is stations 
the rate of 
interfering 
group speed 



will be dominant at the prescribed values 
The values of t and r for which the phase 
ry for a given frequency, therefore, determine 
propagation of this frequency in the mutually 
train of sinusoidal waves, [Ref. 3], called 

c s • 



f (w, r,t ) 



6~>t - k (co) r + 7 T 

n 4 



(2.74) 



= t - r d_ko_ _ o 

d^ dco 



(2.75) 




r _ d co 
T dk„ 



(2.76) 



but 



CO 




Pn 



Therefore 



c = =. 411^) = 

a d k » d k„ 



c i k d c p ^ 

rn n 6 k n (2.77) 



From equation (2.63) we have k n = Ck A / H . Substituting this 
value in equation (2.77), we get 
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(2.78) 



C 




o( d C m 
"do(„ 



From equation (2.71) we can write 



1 



d — c< 


" 1 + (-&.) 1 


" 2 

. 2 


*1 
1 


1 

CL 

X 

1 

X 

o 

i 


d<*„ 2 


l V <V \ 




tf/J 





/ d _ Xo \ 
U dc_ PM = c,XA dc<n # n ; 

V«„ + K 



From equation (2.68) we get 



(2.79) 



d Xn = £ / 1 + S cot X n 

d« n 0 (X„) 



(2.80) 



And then substitution of equation (2.80) and equation (2.79) 
into equation (2.78) yields 





(2.81) 



where 

^x„> = [ 1 + Scot X„ - & X n cos X„/sm 3 X n 1 

(2.82) 
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Figure 2.7 Group and phase speed as function of frequency. 



Figure (2.7) shows the group speed as a function of 
frequency for first three modes. It may be noted that group 
speed is always less than the phase speed. The general 
features for each mode are, [Ref. 2]. 

1 . The first arrival would be nearly sinusoidal waves 
of frequency f n , where f„ denotes the limiting 
frequency for the n th mode. The wave would arrive 
at t = r/C 2 travelling with C e<r> = C 2 . This portion 
of the signal is called ground wave. As the time 
progresses, the frequency and amplitude in the ground 
wave would increase. 

2 . At time t = r/C^ a new train of high frequency 

waves, due to the right hand branch of the group 
velocity curve, would suddenly be superimposed on 
the ground wave travelling with speed Cg n = C., . 

This new high frequency wave is designated as water 
wave. The frequency in the water wave would 

decrease as time progresses while in the ground 
wave the frequency would continue to increase. 

3. Still later f^ and f w would approach each other 
at t=r/Cg 0 , where Cg 0 denote the minimum group speed 
they would coincide. The pressure then would 
consists of a single frequency f A and this portion 
of the signal is called Airy phase. The sequence of 
these events is illustrated in Figure 2.8 . 
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Ill . SOUND PROPAGATION IN A MODELED SHALLOW WATER 

ENVIRONMENT 

A. INTRODUCTION 

The purpose of this experiment was to investigate exper- 
imentally the behaviour of sound propagation in a laboratory 
facility consisting of a layer of water overlaying a layer 
of sand and to compare the results with the normal mode 
theory. 

B. EQUIPMENT 

1. The tank used for experimental work is made of alumi- 
nium plate one centimeter thick. The interior dimensions of 
the tank are : 

Length = 9.95 meters 
Width = 0.94 meters 
Depth = 0.24 meters 

The tank has a carriage which can be used to adjust 
horizontal and vertical positions of receiver. The rails on 
which the carriage moves were leveled to ensure that 
receiver remained at a constant depth as it moved along the 
tank. The leveling arrangement was not precise and it was 
found that the total vertical variation of the receiver 
was within + 1 mm. All the joints of the tank were sealed 
with silicon sealer and the interior of the tank was coated 
with three layers of primer paint (zinc chromate) to avoid 
corrosion. Wooden wedges are fixed on the inner sides of the 
tank to minimize the reflection of sound. The tank was 
filled with water saturated sand to form an 11 cm thick 
bottom layer and a thin layer of water (2 to 3 cm) was added 
above the sand. The sand was degassed by circulating water 
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deep into the sand with a centrifugal pump at moderate pres- 
sure. This process was repeated several times. The top 
surface of the sand was leveled with reference to the water 
surface. This leveled surface of sand was horizonally flat 
within + 1 mm. The layer of water was very shallow and the 
thickness of this layer changes with time due to evaporation 
of the water. To avoid this, the tank was equiped with an 
automatic water level control which maintains the desired 
level of water. 

2. The frequencies used in this work ranged from 25 kHz 
to 100 kHz. The source transducer was built from a rectan- 
gular crystal of piezoelectric material 2.98 cm long and 
1.12 cm wide. The directivity index of the transducer at 50 
kHz is 5.75 dB. The receiver was made from a 4 mm by 4 mm 
piezoceramic cylinder. 

C. MEASUREMENT OF SPEED OF SOUND AND DENSITY 

Figure 3.1 shows the equipment and the circuit used to 
measure the speed of sound in water- saturated sand and in 
water. Two LC-5 transducers were used as source and 
receiver. A high frequency tone burst of the order of one 
MHz was used and the speed of sound was calculated by 
direct measurement of time-of-f light from the oscilloscope 
and the distance between the source and the receiver. 
Experimentally determined values for the speed of sound in 
the sand and in water are shown in Table I and Table II. The 
standard deviation of the measurements is an indication of 
the precision of the measurement (approximately 1 percent). 
The difference between the measured speed of sound in water 
and the value predicted by theory at that temperature is an 
indication of the accuracy of the measurement (approximately 
0.95 percent). The density of the sand was measured by 
weighing 1000 cm 5 of water-saturated sand and calculating its 
mass The measured values are shown in Table I . 
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Figure 3.1 Measurement of speed of sound using 1 MHz tone burst. 



TABLE I 





Speed 


of sound in 


sand 




No. 


Time 


Scale 


Di stance 


C 




Divi sion 


( M. s/div) 


(cm) 


(m/s) 


1 


8.1 


50 


67.8 


1674.07 


2 


6.4 


50 


54.3 


1696.88 


3 


5.5 


50 


45.3 


1647.27 


4 


4.6 


50 


38.0 


1652 . 17 


5 


6.3 


20 


21.1 


1674.6 


6 


8 . 6 


10 


14.2 


1651.16 



Average Speed 


= 


1666.03 


m/s 


Standard Deviation 


- 


17.58 




Percent Deviation 




1.06 


% 



Measurement For Density of Sand 
Mass of the Container = 542 gm 

Mass of 1000 cm 3 of sand + Container = 2650.3 gm 
Density = 2108.3 kg/m 3 
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TABLE II 

O 

Speed of sound in water (Temperature =17 C) 



No. 


Time 

Division 


Scale 
( ^ s/div) 

50 


Di stance 
(cm) 


(m/s) 


1 


3.8 


27.5 


1447.3 


2 


4.5 


50 


33.0 


1466.67 


3 


6.0 


50 


43.5 


1450.0 


4 


7.05 


50 


51.5 


1461.0 


5 


8.8 


50 


64.5 


1465 .91 


6 


5 . 0 


100 


72.2 


1444.0 


7 


5.7 


100 


83 . 7 


1468.42 


8 


6.4 


100 


93.7 


1464.06 



Average Speed of Sound = 

Standard Deviation = 


1458.43 
9 . 1 


m/s 


Speed of Sound from Theory = 


1472.39 


m/s 


Percent Deviation = 


0.95 


% 
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D. VOLTAGE RESPONSE AND SENSITIVITY OF TRANSDUCERS 



The basic characteristics of the source and receiver 
were not known, therefore the voltage response of the trans- 
mitter and the sensitivity of the receiver were measured by 
the reciprocity method (Figure 3.2). The left hand side of 
the Figure indicates the transducers being used as transmit- 
ters and the right hand side as receivers. Experimentally 
measured data are shown in Table III and calculated values 
of voltage response and sensitivity are shown in Table IV. 
The voltage response of the transmitter is plotted in Figure 
3.3 and sensitivity of the receiver in Figure 3.4 . 



45 




46 



TABLE III 



Experimental measurements for 
calibrations of transducers 



F re . 
(kHz) 


V* 

(volt) 


V* 

( volt ) 


Ir 

(amp) 




V* 

( volt ) 


20 


0.00046 


0.0015 


0.0051 




0.00036 


30 


0.00073 


0.002 


0.0075 




0.006 


40 


0.0012 


0.0016 


0.01 




0.012 


50 


0.00175 


0.0074 


0.0127 




0.034 


60 


0.00255 


0.0071 


0.0155 




0.044 


70 


0.003 


0.0227 


0.0185 




0.099 


80 


0.0046 


0.018 


0.0215 




0.224 


90 


0.006 


0.059 


0.025 




0.46 


100 


0.0075 


0.0625 


0.0275 




0.62 


110 


0.01 


0.185 


0.0305 




1 . 15 


120 


0.0123 


0.089 


0.0315 




0.17 


130 


0.016 


0.085 


0.034 




0.193 


140 


0.015 


0.06 


0.037 




0.1 


150 


0.0145 


0.093 


0.043 




0.21 


Di stance 


between source 


and receiver = 


1 


m . 


Applied voltage 




= 


5 


volts 
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TABLE IV 






Voltage response and 


sensitivity 


Fre . 


Sensitivity 


Voltage Response 


(kHz) 


dB ( re V/Pa) 


dB ( re Pa/V @ 1 m) 


20 


-76.37 


5.91 


30 


-68.34 


0.38 


40 


-70.96 


1 . 06 


50 


-63 . 43 


6.84 


60 


-65.78 


8.83 


70 


-59.36 


12.50 


80 


-59.91 


11.03 


90 


-53.95 


15.39 


100 


-54.24 


16.18 


110 


-48.96 


20.32 


120 


-61.86 


26.87 


130 


-63.33 


27 . 94 


140 


-68.10 


29.69 


150 


-63 . 68 


29.07 
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Figure 3.3 Voltage response of transducer. 



o 




jad q.^OA auo ea gp 



n 

X 

X 



>. 

o 

a 

<d 

D 1 

0) 



50 



Figure 3.4 Sensitivity of receiver. 




E. VALIDATION 



To insure that the sound fields set up in the tank would 
be good approximations of the sound fields in the absence of 
sides and bottom (e.g, a shallow water environment of infi- 
nite horizontal extent), measurement were made on the sound 
field set up in the water and in the sand for a source in 
the water layer. 

1 . Sound Pressure as a Function of Receiver Depth 

Figure 3.5 shows the circuit arrangement for the 
measurement of the variation of pressure amplitude with 
range and depth. Continuous waves of the desired frequency 
were generated by a frequency generator and were applied to 
the source through an amplifier. The source was placed at 
mid depth in the water layer which was 2.7 cm thick. The 
signal was received by the receiver and was passed through 
the bandpass filter to eliminate noise. The output from the 
bandpass filter was amplified and displayed on the oscillo- 
scope. To determine the variations of pressure amplitude as 
a function of receiver depth, the transmitter was kept at a 
constant depth in the water layer while the measurements 
were made with the receiver at a fixed range beginning from 
the top surface of the water layer down into the bottom 
layer in steps of 0.5 centimeter. Three frequencies between 
the cut-offs of two first modes were selected for measure- 
ments. The data are tabulated in Tables V, VI, and VII and 
plotted in Figures 3.6, 3.7, and 3.8 . The experimental and 
theoretical results agree satisfactorily within the water 
layer. Near the cut-off frequency of the propagating mode, 
the experimental amplitude appears to decay more rapidly 
with depth in the bottom layer than theory predict. Decay is 
slower for intermediate and high frequencies. 
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Figure 3.5 Measurement of pressure variations with range and depth. 




TABLE V 



Sound pressure as function of depth: 
f = 35 kHz, n = 1, L = 1.5 m 



No. 


Depth( cm) 


Volts 


Normal 


1 


0.0 


0.012 


0.12 


2 


0.5 


0.028 


0.28 


3 


1.0 


0.056 


0.56 


4 


1.5 


0.081 


0.81 


5 


2.0 


0.097 


0 . 97 


6 


2.5 


0. 1 


1.0 


7 


3 . 0 


0.097 


0.97 


8 


3.5 


0.083 


0.83 


9 


4.0 


0.067 


0.67 


10 


4. 5 


0.052 


0.52 


11 


5.0 


0.042 


0.42 


12 


5.5 


0.041 


0.41 


13 


6.0 


0.038 


0.38 


14 


6.5 


0.034 


0.34 


15 


7.0 


0.028 


0.28 


16 


7.5 


0.023 


0.23 


17 


8.0 


0.018 


0. 18 


18 


8.5 


0.015 


0.15 


19 


9 . 0 


0.016 


0. 16 
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Depth (cm. 



Normalize Amplitude 



0.0 0.2 0.4 0.6 0.8 1.0 




Figure 3.6 Pressure variation as function of receiver 

depth: f = 35 kHz n = 1. 
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No 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 



TABLE VI 



pressure as 
= 50 kHz, n 

Depth( cm) 
0.0 
0.5 
1.0 

1.5 
2.0 

2 . 5 

3.0 

3 . 5 

4.0 

4.5 

5.0 

5.5 

6.0 

6.5 

7.0 

7.5 

8.0 

8.5 
9 . 0 



function of 
= 1, L = 1.5 

Volts 

0.04 

0.21 

0.41 

0.58 

0.67 

0.65 

0.58 

0.435 

0.34 

0.245 

0.15 

0.1 

0.063 

0.05 

0.05 

0.047 

0.027 

0.02 

0.04 



depth : 
m 

Normalized 

0.06 

0.31 

0.61 

0.87 

1.0 

0.97 

0.87 

0.65 

0.51 

0.37 

0.22 

0.15 

0.09 

0.07 

0.07 

0.07 

0.04 

0.03 

0.06 
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Normalize Amplitude 



0.0 0.2 0.4 0.6 0.8 1.0 




Figure 3.7 Pressure variation as function of receiver 

depth: f = 50 kHz n = 1. 
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TABLE VII 





Sound 

f 


pressure as 
= 70 kHz, n 


function 
= 1, L = 


of depth: 
2.0 m 


No. 




Depth( cm) 


Volts 


Normalized 


1 




0.0 


0.095 


0.05 


2 




0.5 


0.53 


0.28 


3 




1.0 


1.15 


0.6 


4 




1.5 


1.7 


0.89 


5 




2.0 


1 . 9 


1.0 


6 




2 . 5 


1.75 


0.92 


7 




3.0 


1.2 


0 . 63 


8 




3.5 


0.7 


0.37 


9 




4.0 


0.24 


0.13 


10 




4.5 


0.205 


0.11 


11 




5.0 


0.2 


0.1 


12 




5.5 


0.195 


' 0. 1 


13 




6.0 


0.13 


0.07 


14 




6.5 


0.043 


0.02 


15 




7.0 


0.06 


0.03 


16 




7.5 


0.065 


0.03 


17 




8.0 


0.1 


0.05 


18 




8.5 


0.105 


0.055 


19 




9.0 


0.073 


0.04 
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Normalize Amplitude 

0.0 0.2 0.4 0.6 0.8 1.0 




Figure 3.8 



Pressure variation as 
depth: f = 70 kHz n 



function of 
= 1 . 



receiver 
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2 . Sound Pressure as a Function of Range 

To measure the variation of peak pressure amplitude 
with range, the transmitter was kept at constant depth and 
the receiver's depth was adjusted slightly at each range to 
measure peak amplitude at that range. Precautions had to be 
taken to ensure that the water surface is disturbed as 
little as possible because variation in water depth had a 
significant effect on the pressure amplitude. Measurements 
were made for three frequencies between the cut-offs of 
the first two modes. The data are tabulated in Tables VIII, 
IX, and X and is plotted in Figures 3.9, 3.10, and 3.11 
respectively. Assuming a simple model for which the loss in 
the bottom effects that portion of the energy density of the 
normal mode found in the bottom, which then "robs" the 
energy density in the water layer, then the effective 
absorption coefficient can be estimated as (Ref. personal 
communication with Prof. Alan B. Coppens) 



Where Ab/A is fraction of the energy density of the normal 
mode in the bottom and 1/4F (F is the frequency in kHz) 
yields the expected attenuation to be found in the sand used 
for the bottom [Ref. 4]. Comparison of the measured data 
with the simplestic model is very good, suggesting that 
indeed the sound field is subjected to minimal effects from 
the side walls. 



a 



_L FA* 
4 A 
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TABLE VIII 



Sound pressure as function of range: 
H = 2 cm f = 50 kHz n = 1 



No. 


Distance(m) 


Volts 


Normalized 


1 


1 . 00 


0.45 


1.0 


2 


1.25 


0.34 


0.756 


3 


1.50 


0.255 


0.567 


4 


1 . 75 


0.195 


0.433 


5 


2 . 00 


0.2 


0.444 


6 


2.25 


0. 168 


0.373 


7 


2.50 


0. 108 


0.24 


8 


2.75 


0.155 


0.344 


9 


3.00 


0. 105 


0.233 


10 


3 . 25 


0.078 


0.173 


11 


3.50 


0.022 


0.049 


12 


3.75 


0.047 


0.104 


13 


4.00 


0.059 


0.131 


14 


4.25 


0.057 


0.127 


15 


4.50 


0.043 


0.095 


16 


4.75 


0.042 


0.093 


17 


5.00 


0.031 


0.069 
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Figure 3.9 Sound pressure variation as function of range: f = 50 kHz 



TABLE IX 



Sound pressure as function of range: 
H = 2 cm f = 70 kHz n = 1 



No. 


Distance(m) 


Volts 


Normalized 


1 


1.00 


2.08 


1.0 


2 


1.25 


1 . 6 


0.769 


3 


1.50 


1.7 


0.817 


4 


1.75 


1.0 


0.481 


5 


2.00 


0.93 


0.447 


6 


2 . 25 


1.1 


0.529 


7 


2.50 


0.85 


0.409 


8 


2 . 75 


0.68 


0.327 


9 


3.00 


0.87 


0.418 


10 


3.25 


0.62 


0.298 


11 


3.50 


0.17 


0.082 


12 


3 . 75 


0.61 


0.293 


13 


4.00 


0.58 


0.279 


14 


4.25 


0.36 


0.173 


15 


4.50 


0.315 


0.151 


16 


4.75 


0.3 


0.144 


17 


5.00 


0.29 


0.139 
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Figure 3.10 Sound pressure variation as function of range: f = 70 kHz 



TABLE X 



Sound pressure as function of range: 
H = 2 cm f = 100 kHz n = 1 



No. 


Distance(m) 


Volts 


Normalize 


1 


1.00 


1.95 


1.0 


2 


1 . 25 


1.4 


0.718 


3 


1.50 


1.4 


0.718 


4 


1.75 


1 . 05 


0.538 


5 


2 . 00 


1.0 


0.513 


6 


2.25 


0.88 


0.451 


7 


2.50 


0.83 


0.426 


8 


2.75 


0.85 


0 . 436 


9 


3.00 


0.74 


0.379 


10 


3.25 


0.63 


0.323 


11 


3 . 50 


0.7 


0.359 


12 


3 . 75 


0.28 


0.144 


13 


4.00 


0.65 


0.333 


14 


4.25 


0.5 


0.256 


15 


4.50 


0.41 


0.21 


16 


4.75 


0.49 


0.251 


17 


5 . 00 


0.41 


0.21 
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Figure 3.11 Sound pressure variation as function of range: f = 100 kHz 



Sound Pressure Variation Across the Width of the 



3 . 

Tank 

Sound pressure variations were measured across the 
width of the tank at a one meter distance from the source 
for three different frequencies. The measured data shown in 
Table XI, XII, and XIII and plotted in Figures 3.12, 3.13, 
and 3.14. The results show slight effects of interference 
either from the walls or from the bottom of the sand layer. 
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TABLE XI 



■ansverse 


measurement 
f = 50 kHz, 


of sound pressure: 
L = 1 m 


. stance 
(cm) 


Volts 


Distance 

(cm) 


Volts 


-32 


0 . 38 


-3 . 5 


0.4 


-30.5 


0.215 


-1.7 


0 . 55 


-29 


0 . 38 


0 


0.40 


-27.5 


0.34 


2 


0 . 55 


-26.4 


0.39 


4 


0.4 


-25.2 


0.32 


6.7 


0.53 


-23 . 7 


0.47 


8.5 


0.34 


-22 


0.33 


10.7 


0.49 


-20 


0.5 


13 . 8 


0.36 


-18 


0.33 


16.3 


0 . 49 


-16.3 


0.5 


18.8 


0.32 


-15 


0 . 42 


20.2 


0.4 


-13 . 5 


0.5 


22.2 


0.31 


-11 . 8 


0.36 


23 . 5 


0.38 


-9.9 


0.53 


24.7 


0.33 


-7.2 


0 . 42 


26 


0.37 


-5.3 


0.53 


27 . 7 


0.25 
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Figure 3.12 Transverse variation of sound pressure: f = 50 kHz, 




TABLE XII 



Transverse 


measurement 
f = 70 kHz, 


of sound pressure: 
L = 1 m 


Distance 

(cm) 


Volts 


Di stance 
(cm) 


Volts 


-32 


1.7 


-1.3 


1.7 


-30.6 


1 . 25 


0.4 


2.47 


-29.4 


1.4 


2 


1.6 


-28.6 


1.25 


3 . 7 


2.17 


- 27.2 


1.8 


5.2 


1 . 35 


-26 


1.25 


7 


2.23 


-24.9 


1 . 75 


8.7 


1 . 45 


-23 . 7 


1.2 


10 


2 . 12 


-22 . 5 


2.05 


11 


1 . 85 


-21.3 


1.2 


12.3 


2.2 


-19 . 8 


2.2 


13 . 7 


1 . 4 


-18.4 


1.5 


15 


2 . 1 


-16 . 8 


1.9 


16.5 


1 . 65 


-15.6 


1.7 


17.8 


1.9 


-14.2 


2 . 1 


18.8 


1.57 


-13.3 


1.82 


19.8 


1.75 


-12 


1.98 


21.2 


1.3 


-10 . 7 


1.75 


22.5 


1.62 


-9 . 2 


2 . 2 


23 . 6 


1 . 35 


-7.8 


1.55 


24.5 


1.6 


-6 


2 . 5 


26 


1 . 0 


-4.5 


1 . 65 


27 . 8 


1.5 


-3 


2 . 32 
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Figure 3.13 Transverse variation of sound pressure: f = 70 kHz, 



TABLE XIII 



Transverse 


measurement 
f= 100 kHz, 


of sound pressure: 
L = 1 m 


Distance 

(cm) 


Volts 


Di stance 
(cm) 


Volts 


-32 


1.25 


-4 


1.47 


-30.8 


0.82 


-1.9 


1.9 


-29.6 


1.5 


-0.4 


1 . 52 


-28.4 


0.94 


0.5 


1.7 


-27 


1 . 4 


1 . 5 


1.6 


-25.5 


0.94 


2.4 


1.7 


-23.8 


1 . 62 


3 . 4 


1.62 


-22.6 


1.15 


5.3 


1.22 


-20.5 


1 . 76 


7.8 


1 . 67 


-18.7 


1 . 55 


10 


1 . 1 


-17.5 


1.8 


12 . 5 


1.35 


-16.4 


1 . 5 


15 . 6 


0.82 


-15.4 


1.9 


16.5 


1.22 


-14.4 


1 . 5 


18 


0.63 


-13 . 5 


1.8 


19.8 


1.05 


-13 


1 . 6 


21.2 


0.67 


-11 . 7 


2 . 2 


22.5 


0.95 


-10.7 


1.65 


24.7 


0.49 


-9.5 


1.8 


26 


0.68 


-8.7 


1 . 65 


27 . 5 


0.42 


-7 


1.88 







71 



o 




( 9pnq.ixdiuv * jon ) Box 



£ 

rH 

II 



72 



Figure 3.14 Transverse variation of sound pressure: f = 100 kHz, 




F. SUMMARY AND CONCLUSION 



In the preceding sections, the experimental curves for 
the pressure amplitude as a function of depth and range for 
three different frequencies for which only the first normal 
mode propagates have been compared with normal mode theory 
and a simple model for attenuation of the normal mode. The 
depth dependent curves both in the water and in the sand 
agrees satisfactorily with the theory. The range dependence 
obtained in the water exhibit excellent agreement with the 
simple modeling of losses resulting from the relatively 
highly attenuating bottom. Measurements both in range and 
transverse dimensions suggests that the sound field is 
subjected to minimal effects from the side walls. The 
leveling arrangement needs improvement 
frequencies minor alterations in water 
pressure amplitude considerably. 



because at higher 
depth changes the 
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